Effect of Dilution on Spinodals and Pseudospinodals 
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We investigate the effect of quenched dilution on the critical and spinodal points in the infinite 
range (mean-field) and long-range (near-mean-field) Ising model. We find that unlike the short-range 
Ising model, the effect of the dilution is not simply related to the divergence of the specific heat, i.e., 
the Harris criterion. We also find the mean-field behavior differs from that of d = 4 (upper critical 
dimension) nearest neighbor model at the critical point. These results are an important first step 
for understanding the effect of the spinodal on the nucleation process as well as the properties of the 
metastable state in systems of considerable interest in material science, geophysics, and econophysics 
which in general have defects. 
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Long range interactions are common as well as impor- 
tant in both physical and social systems. Metals 0, Q 
and earthquake faults [3(, both with long-range elastic 
interactions, as well as biological and polymer materi- 
als 0] all have processes strongly influenced by the long 
range nature of their interactions. In addition, economic 
systems as well as societal interactions that facilitate 
the transmission of diseases [H[ have effective long range 
"interactions" due to globalization. In all of these sys- 
tems metastable states are common [f| 0| and nucleation 
plays an important role. In addition, it has been estab- 
lished that in long range systems the nucleation process 
is strongly affected by the pseudospinodal @- It is im- 
portant to distinguish mean-field (mf) systems with in- 
finite range interactions 0, [l(| from systems with long 
but finite range interactions. We will refer to the latter 
as near- mean- field (nmf). Mean- field systems have spin- 
odals while nmf systems have pseudospinodals jia.lll, 12 j. 
Pseudospinodals behave as spinodals unless they are ap- 
proached too closely as we will discuss below. 

Although the spinodal and pseudospinodal have a 
strong effect on the metastable state and nucleation in 
defect free mf and nmf systems [4, 8], the systems men- 
tioned above generally have defects. In materials such 
as metals, the defects can take the form of impurities 
and/or vacancies In economic and other social sys- 
tems agents that behave differently than the majority 
can be thought of as defects. Examples include the pres- 
ence of immune individuals in models of disease trans- 
mission Q and agents that do not participate in wealth 
exchange in economic models. This raises an important, 
but as yet unanswered, question: How, if it all, does the 
presence of defects affect the spinodal and pseudospin- 
odal? This question is analogous to what was asked 
about critical points in systems with defects and short 
range interactions several decades ago. In the 70's, Har- 
ris 13[ developed a criterion associated with the specific 



short range Ising models in the presence of vacancies. In 
this paper we will investigate the same question for mf 
and nmf critical points, spinodals and pseudospinodals. 
As we will show, mf and nmf Ising systems have a more 
complicated response to defects than nearest neighbor 
models. Moreover, the response is different than what 
one finds in d — d c =4, the Ising upper critical dimen- 
sionBlll. 

We begin by noting that the Ginzburg criterion for mf 
to be a good approximation in spin systems is that the 
fluctuations in the magnetization are small compared to 
its mean value. At the critical point Q 



G~ 



t d XT 



->■ G = R d e 



,d 2-d/2 



— > OO, 



(1) 



where for interactions with range R, the correlation 
length is £ = Rer 1 / 2 , the reduced temperature is e = 
(T — T c )/T c , and T c is the critical temperature. The sus- 
ceptibility, \T oc e _1 , and the magnetization per spin, 
4> oc e 1 / 2 . We can use G to distinguish between mf 
(G — > oo) and nmf (G> 1 but not infinite) for all d. 

Near spinodals or pseudospinodals the same consider- 
ations lead to the criterion 81 



G s = R d Ah 3 / 2 ~ d / 4 » 1 



(2) 



heat exponent in the defect free, or pure, system that pre- 
dicted whether the critical point remaines well defined in 



for there to be a well pronounced pseudospinodal; G s — > 
oo for a true spinodal [1, [H| . Here the correlation length 
is £ = RAh^ 1 / 4 , the susceptibility is \t ~ A/i -1 / 2 , and 
the difference between the magnetization and its value 
at the spinodal is A<f> ~ Ah 1 / 2 where Ah = (h s — h)/h s 
where h s is the value of the magnetic field at the spinodal. 
Some of the simulations we discuss below involve fully 
connected models, i.e., every spin interacts with every 
other spin. As in finite size scaling, we take the volume 
of the system, TV, to be either £ d = R d t~ d / 2 (critical 
point) or R d Ah- d ' i (spinodal). Then G = Ne 2 (critical 
point) and G s — N Aft, 3 / 2 (spinodal). 

The Harris criterion [l3| (HC) compares the width of 
the distribution of vacant sites (p = fraction of occupied 
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sites) to the distance from the continuous transition. At 
the critical point this leads to no change in exponents if 
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Using again the fact that £ 
G = R d e 2 - d / 2 » 1. 

How does the HC relate to the specific heat in mf and 
nmf systems? We can obtain the specific heat by calcu- 
lating the energy fluctuations per unit volume. In mf and 
nmf, the energy fluctuations are the square of the mag- 
netization fluctuations. The magnetization fluctuations 
per unit volume in mf and nmf do not scale as e" (with 
P = 1/2 at mean- field critical points) but rather as [|| 
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At critical points, the scaling in Eq. (0} can be seen 
by looking at a Landau-Ginzburg- Wilson Hamiltonian 
(LGW) in zero field 15| We scale all lengths with £ and 
m(x) by e 1 / 2 and obtain 
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where y is the scaled length and 4>(y) is the scaled mag- 
netization. As can be readily seen from using Eq. ([5]) 
to calculate the partition function, 4>(y) is of order 
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This scaling reduces the LGW Hamil- 



tonian to a Gaussian because we can ignore the scaled 
<f) A term. For T > T c and h — 0, the specific heat at 
constant h is 



C\ 



,1/2 



y/ R d e 2-d/2 



R d e -d/2 



1 



R d £ 2-d/2 



(G) 



This is the square of the energy fluctuations per unit vol- 
ume and hence for T > T c , the HC for the mf (or nmf) 
critical point to be unaffected by the dilution becomes 
Cv <C 1. In Fig. Q] we plot the energy fluctuations (spe- 
cific heat) of the non-diluted model as a function of e in 
d = 1, 2 (blue and red, respectively). For T > T c the di- 
vergence is consistent with Eq. ([5]), namely e~ 3 / 2 in d = 1 
and (T l in d = 2. 

There are several points to be made. The first is that 
at the critical point the nmf criterion that G 3> 1 is the 
same as the HC so that if the system is well approximated 
by mf theory, the critical point remains well defined and 
unchanged in the presence of vacancies. In addition, for 
G — > oo (mf) the critical point is unchanged by the dilu- 
tion. The second point is that for R 3> 1 but finite, the 
affect of the dilution depends on the distance from the 
critical point e. Hence dilution in the d = 3 Ising model 
has no affect on the apparent mean-field critical point 
if G > 1 and finite. However, as e decreases and the 
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FIG. 1: (color online) Plots of the fluctuations in the en- 
ergy(specific heat) as a function of the distance from the crit- 
ical point for the pure system in d = 1 and 2 for T > T c . 
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FIG. 2: (color online) Susceptibility and specific heat as a 
function of e for a system with average dilution of 10%. Sites 
were designated as empty with a 0.1 probability. Twenty 
systems were simulated. 



system crosses over to d = 3 short-range Ising behavior, 
the vacancies alter the critical behavior from what one 
obtains in the pure model. 

Our result is somewhat different than one obtained by 
Ballesteros et al. [13] who looked at the nearest neigh- 
bor dilute Ising model in d = d c = 4 where the specific 
heat exponent a = 0, so the HC is indeterminate [l3j ]. 
Despite no power-law divergence, the pure system does 
have a logarithmically diverging specific heat. It was 
found in Ref. [14J that the logarithmic corrections to the 
exponents changed; however, the exponents themselves 
did not. For the mf (G — > oo) system we consider, the 
HC predicts that the nature of the critical point remains 
the same as the pure system. 

In Fig. [3] we plot the susceptibility \ and specific heat 
Cy as a function of e > for the system with dilution and 
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G ^> 1 but finite for the fully connected model. As we 
approach the critical point and G decreases, the spread 
in both Cy and \ increases consistent with the HC. (See 
Fig.EJ) 

We next turn to the spinodal. The HC is 



[£ rf p(l — p)] 1 / 2 constant 

With f = RAh- 1 / 4 , the HC becomes 
Ah- 1 / 2 



< Ah. 
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The specific heat is again related to the fluctuations of 
the energy. Near the spinodal the magnetization is non- 
zero. By using essentially the same argument as outlined 
in Eq. ([5|) , the fluctuations in the magnetization scale as 
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where we have expanded the LGW Hamiltonian around 
the magnetization at the spinodal 0. 

The specific heat (energy fluctuations) at the spinodal 

is 
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Note that near the spinodal the magnetization does not 
go to zero as the spinodal is approached. Near the pseu- 
dospinodal, the HC is not simply related to the specific 
heat. What is required for the pseudospinodal to be un- 
changed by the presence of dilution according to the HC 
is that the specific heat divided by the Ginzburg parame- 
ter be small. This result implies that in mean-field where 
the R — > oo before Ah — > 0, [l(| the spinodal is unaf- 
fected by the dilution. In Fig.[3]we plot the specific heat 
(energy fluctuations) in the long-range Ising model with- 
out dilution as the spinodal is approached. The result is 
consistent with Csp ~ Ah^ 1 / 2 . 

In Fig. 3] we plot 20 realizations of the isothermal sus- 
ceptibility and specific heat as a function of h near the 
pseudospinodal for the fully connected model with 9 x 10 4 
spins. The empty sites are distributed with a probability 
of 0.10 as for the critical point. Note that the curves for 
each simulation overlap until Ah becomes too small at 
which point the pseudospinodal is smeared as evidenced 
by the fact that the curves show a distinct spread. This is 
consistent with our analysis. For both the critical point 
and the pseudospinodal the spread caused by the dilution 
in the fully connected model narrows for a given value of 
the distance from the pseudospinodal or the critical point 
as the system size is increased and hence R increases as 
would be expected from our analysis. 
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FIG. 3: The specific heat for the pure system plotted as a 
function of h for d = 2. The fit is consistent with a divergence 
with exponent 1/2 consistent with the scaling analysis. 
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FIG. 4: (color online) The susceptibility and specific heat as 
a function of the distance from the spinodal for the dilute 
fully connected model. Twenty systems with a probability 
of a site being diluted of 0.1 were simulated. The existence 
of the spread for small Ah is consistent with the theoretical 
analysis. 



In Fig. [5] we plot the width or spread of the suscepti- 
bility near the critical point (normalized by the suscepti- 
bility) for a 20% and 40% dilution in the fully connected 
model for two system sizes. The width increases as e de- 
creases as expected. We did not display data for e < 0.04 
because the nmf picture breaks down for e ~ 0.05. Simi- 
lar results are obtained for the pseudospinodal 
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cause the interaction range is a function of the system 
size, we expect the spread of the susceptibility associated 
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on nucleation and pseudospinodals in more complicated 
models associated with materials such as the long-range 
Ising antiferromagnet Q? What is the effect of other 
types of damage such as fixed spins or dislocations? 
What is the effect of various forms of damage and impu- 
rities in models in other fields such as econophysics and 
the spread of disease? The answers to these questions are 
of considerable importance in a variety of fields of study 
ranging from materials and biophysics to economics and 
social science. 

This research was supported by the Department of En- 
ergy through grant DE-FG02-95ER14498. We acknowl- 
edge useful conversations with H. Gould, J. Silva and J. 
B. Rundle. 



FIG. 5: (color online) Width of the susceptibility (variance 
normalized by the mean) distribution near the critical point 
(see Fig. HJ for (a) q = 0.20 and (b) q = 0.40 with a cutoff of 
e = 0.04. 



with the dilution to decrease as the system size increases 
for the same value of e. This is consistent with the scaling 
analysis. A similar result is found near the spinodal [l7j . 

In conclusion we have found that the effect of dilution 
in mf and nmf systems on critical points is more subtle 
than the effect in systems with short range interactions 
below the upper critical dimension, and it is different 
from the behavior of short range systems in d — d c di- 
mensions. We also found that in mf systems the spinodal 
is not changed by dilution but in more physically realistic 
nmf systems the spread and consequent rounding of the 
pseudospinodal depends on the distance from the point 
of the apparent divergence. At the spinodal the HC is 
not simply related to the mf specific heat exponent [l6j |. 
Near T c for h = the Harris criterion is simply related 
to the specific heat for T > T c . Near the critical point 
but for T < T c we found a result similar to what one 
finds near the spinodal. These results will be published 
in Ref. (T^j. 

Finally, our results indicate that the effect of the spin- 
odal on nucleation in a wide variety of systems is not 
eliminated by the presence of dilution and that spinodal 
nucleation remains an important feature of systems with 
long range interactions. 

Our results raise many questions and suggest several 
new research directions. What is the effect of dilution 
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